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Experimental data from the Mega Amp Spherical Tokamak (MAST) is used to show that the 
inverse gradient scale length of the ion temperature RjhT i (normalized to the major radius R) has 
its strongest local correlation with the rotational shear and the pitch angle of the magnetic field (or, 
equivalently, an inverse correlation with q/e, the safety factor/the inverse aspect ratio). Furthermore, 
RjhTi I s found to be inversely correlated with the gyro-Bohm- normalized local turbulent heat flux 
estimated from the density fluctuation level measured using a 2D Beam Emission Spectroscopy 
(BES) diagnostic. These results can be explained in terms of the conjecture that the turbulent system 
adjusts to keep R/LT t close to a certain critical value (marginal for the excitation of turbulence) 
determined by local equilibrium parameters (although not necessarily by linear stability). 



Introduction. A key physics challenge posed by mag- 
netically confined plasmas in fusion devices is how their 
internal energy can be kept from being transported too 
fast from the core to the periphery. The problem is pri- 
marily one of turbulent transport, the temperature gra- 
dient between the edge and the core of a toroidal plasma 
supplying the free energy for the turbulent fluctuations 
that then enhance the effective thermal diffusivity and 
relax the gradient. It is the ion temperature gradient 
(ITG) that is expected to cause the most virulent insta- 
bilities (on ion Larmor scales) and then to be self- 
consistently limited by the resulting turbulence Q • If we 
view the edge ion temperature as fixed by the physics and 
engineering aspects of the tokamak design that will not 
concern us here the key question is how to maximize 
the ion temperature gradient. We therefore wish to in- 
quire, experimentally, what this gradient depends on and 
how. Motivated by the fact (or the conjecture) that the 
ion-scale turbulence is largely determined by the local (to 
a given flux surface) equilibrium conditions and in 
turn acts back to adjust them locally, we ask what local 
parameters are most strongly correlated with the corre- 
sponding value of R/LTi, the inverse radial gradient scale 
length of the ion temperature (L^,. = \d1nTi/dr\) nor- 
malized to the major radius R of the torus. 

How universal are any such measured dependences 
likely to be for situations with different global condi- 
tions, e.g., different neutral-beam-injection (NBI) heat- 
ing powers? It has been recognized for some time 
that the turbulent heat flux tends to increase very 
strongly (much faster than linearly) with R/L^, a phe- 



nomenon known as "stiff" transport [llHla ]. If (or 
when) the transport is stiff, any experimentally mea- 
sured relationship between R/L^ and other equilibrium 
parameters should be quite close to some critical man- 
ifold in the parameter space separating dominant tur- 



bulent transport from a non-turbulent or weakly turbu- 
lent state (the "zero-turbulence manifold" [HI). This 
critical manifold would be independent of the power in- 
put and can be represented as a local parameter depen- 
dence of the critical temperature gradient R/Lr i C = 
f(q,s,s,U' <j> ,R/L n ,R/L TB ,i/ ii ,Ti/T e ,j3i,...), where q is 
the safety factor (number of toroidal revolutions per one 
poloidal revolution of the magnetic field around the torus 
on a given flux surface), e = r/R the inverse aspect ratio 
(r is the minor radius of the flux surface) , s = dlnq/dhir 
the magnetic shear, UL = dU^/dr the radial shear of the 
mean toroidal rotation velocity U^, L n and Lt c the gra- 
dient scale lengths of the plasma density and electron 
temperature, v>u the ion collision rate, Ti/T e the ion-to- 
electron temperature ratio, /3, = %imTi/B 2 the ion-to- 
magnetic pressure ratio and "..." stand for everything 
else (e.g., the many parameters required to fully describe 
the magnetic configuration) [ijj. We stress that R/L^.c 
need not be the same as the threshold for the existence 
of linearly unstable eigenmodes. Two known examples 
when it is not are the "Dimits upshift" of R/L^ above 
the linear stability threshold |18l - t20l ] and the case of suf- 
ficiently large when the system is linearly stable but 
strong transient excitations [21 , 22 1 lead to sustained sub- 
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-251]. Thus, in general, there is a 
nonlinear threshold with some definite dependence on lo- 
cal equilibrium parameters. 

Recent theoretical [21, 22 1 and numerical [lj| investiga- 
tions suggest that q/e and may be the most important 
such parameters, at least at low s. Let us explain why 
this is. It i s well known, both from experimental measure- 
13 and theory [1, [H, [IH , that strong (finite- 



ments 



Mach) flows in a tokamak are predominantly toroidal 
(certainly when plasma is heated by tangential neutral 
beams, which produce a toroidal torque) [301 ] - Therefore, 
any radial shear in the toroidal flow results in sheared 
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flow in both the perpendicular (U'j_ = (B p /B) UL, B p is 
the poloidal field) and parallel (f7,j = (B^/B) U^, B^ is 
the toroidal field) directions. While perpendicular flow 
shear is known (theoretically [23, 24, 31-37} and experi- 
mentally fili - fljl . Islj - Eol ] ) to suppress turbulence and the 
associated transport, parallel flow shear can drive turbu- 
lence via the ''parallel-velocity-gradient" (PVG) instabil- 
ity 



21, 22, 41 



The average ratio of these two shearing 
rates on a flux surface, U'^/U'^ = B^/B v , can be approx- 
imated by q/e and so the degree to which sheared equi- 
librium flow suppresses or drives turbulence is expected 
to depend on this parameter. Indeed, numerical studies 
of ITG- and PVG-driven turbulence have shown that the 
critical R/L^.c at any given U'^ increases with decreas- 
ing q/e (at least for low s [IH); while at any given q/e, 
R/LTi.c increases with increasing U'^ provided the latter 
is not too large 111 23-25| (as in most real tokamaks). 

A comprehensive numerical parameter scan of the de- 
pendence of R/ Lx it c on all other potentially important 
local quantities (s, R/L n , R/Lt b , Vu, Ti/T e , Pi, etc.) is 
probably unaffordable in the near future. Faster progress 
can be made experimentally. In this Letter, our first 
goal is to establish, based on a relatively sizable dataset 
for MAST, what the most important parameters for the 
critical manifold arc: we will show that, indeed, the lo- 
cal value of R/L^ is most strongly correlated inversely 
with the local q/e and positively with the local rotational 
shear — consistently with the result obtained in (lij . 

Our second goal is to obtain an experimental signature 
that the measured R/ is determined by — or, more 
precisely, correlated with — the local characteristics of 
the ion-scale turbulence, directly measured by the 2D 
beam emission spectroscopy (BES) diagnostic [42j. We 
will show that not only does a strong correlation between 
R/L^ and an estimated turbulent heat flux level exist 
but its (at the first glance, counterintuitive) inverse na- 
ture is consistent with R/Lt ( staying close to the critical 
threshold R/L^.c and hence with stiff transport. 

Equilibrium parameters. A database was compiled of 
equilibrium quantities (and turbulence characteristics; 
see below) from 39 neutral-beam-heated discharges from 
the 2011 MAST experimental campaign. These dis- 
charges had a double-null diverted (DND) magnetic con- 
figuration, no pellet injection and no applied resonant 
magnetic perturbations. Mean electron density n e and 
temperature T P were measured with the Thomson scat- 
tering system [43j, mean impurity ion (C 6+ ) temperature 
Ti and the toroidal flow velocity with the Charge eX- 
change Recombination Spectroscopy (CXRS) system [44 1 
(we assumed that in these discharges the impurity and 
bulk ions have negligible differences in their temperature 
and flow velocities [45j). The local magnetic pitch angle 
{B p /B$) was measured with the Motional Stark Effect 
(MSE) diagnostic [46|; pressure- and MSE-constrained 
EFIT equilibria [47| were used to obtain the field strength 



B. All parameters were determined over 5 ms intervals 
either by averaging if the diagnostic's temporal resolu- 
tion was smaller or by interpolation if it was larger than 
5 ms. Only data points from a limited range of minor 
radii 0.6 < r/a < 0.7 (r = a is the edge of the plasma) 
were used, in order to minimize any correlations between 
various quantities due to their profile dependence alone 
(thus, we did not attempt to prove locality here; see, how- 
ever, In total, 988 data points were available (out 
of which 109 points are from H-mode discharges). 

From this information, we constructed 7 local di- 
mensionless parameters, which, motivated by theoretical 
models and common sense, we deemed a priori the most 
important ones (we also give the range of variation of 
each parameter): R/L T% € [0.08,20.3], q/e G [4.0,16.3], 
s e [1.2,6.0], 7 £ = U' ± T Bt = nrU'Jv thi € [0.005,2.5], 
R/L n e [0.04,13.8], R/L Tf e [1.43,22.7], v« = v u r at e 
[0.003,0.12], T l /T e e [0.5, 1.7]. The ion collision rate v u 
and the perpendicular velocity shear U' , (which is used 
instead of U'^) were normalized, as in [id ], to the ion par- 
allel streaming time r st = A/vthi, where vthi = y/2Ti/rrii 
and A = -KrB/B v is the connection length (the approx- 
imate distance along the field line from the outboard to 
the inboard side of the torus, expected to determine the 
parallel correlation scale of the turbulence [l2| ; if the 
flux surfaces had been circular, A nqR). The local 
magnetic configuration is represented by q/e and s. The 
choice of q/e was motivated by the physical considera- 
tions outlined in the Introduction; since e varied little 
in our database, we cannot distinguish any individual 
correlations of R/L^ with q and e. It is left for fur- 
ther study whether other properties of the flux surfaces 
matter (e.g., Shafranov shift, triangularity, elongation, 
etc.; some of these may, in fact, affect the stiff-transport 
threshold El]). We have not included n, T h T e , 



which arc not normalizablc by any natural local quan- 
tities; note that R/L^ usually has a large but trivial 
correlation with Tf. larger temperature gradients lead to 
larger temperatures in the core. We also have excluded 
Pi = %imTi/B 2 because, in the absence of large variation 
of B in our dataset, Pi is simply the normalized ion pres- 
sure and, similarly to Ti, has a large positive correlation 
with R/L^ (it remains to be investigated whether larger 
magnetic fluctuations at larger pi are large enough to 
have a nontrivial effect on turbulent transport 491-551]). 

Correlation analysis. We perform a Canonical Corre- 
lation Analysis (CCA) [HI] with R/L^ treated as the de- 
pendent variable and the other 7 local parameters item- 
ized above as independent ones. This amounts to finding 
the maximum correlations between \rt{R/ L^) and linear 
combinations of logarithms of 1, 2, 3, . . . , or 7 other pa- 
rameters, leading to an effective statistical dependence 
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TABLE I. Results of CCA performed assuming Eq. {TJ. 
Wherever appears, that means that the CCA was performed 
without including the corresponding parameter. 
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This is of course not valid if the dependence of R/L^ 
on any of the parameters is non-monotonic. A non- 
monotonic dependence on 7^ is, in fact, expected, with 
R/L^ first increasing, then decreasing at larger values of 
je due to increased transport from the PVG-driven tur- 
bulence [l6, 23 - 25 1 . However, the range of values of je 



in our database do not extend to sufficiently high values 
for such a dependence to be observed (see Fig.Q]). 

The results are shown in Table HI where the values of 
the canonical correlation (i.e., the correlation coefficient 
between the logarithms of R/L^ and the right-hand of 
Eq. are given together with the corresponding expo- 
nents ai, . . . , 017. We start by calculating the individual 
correlations of R/ with each of the 7 parameters and 
then include pairs, triplets, etc., only if the correlation 
improves. We see that the strongest individual correla- 
tion of R/L Ti arc with q/e (61%) and j E (46%). The 
overall fit is measurably improved (66%) if both are in- 
cluded. Including further parameters does not make a 
significant difference; the third strongest (although not 
very strong) dependence is on 

The dependence of R/L^ on q/e and te is shown in 
Fig.Q] R/ gen erally increases with decreasing q/e and 



increasing 7^ 57] . broadly consistent with the expecta- 



tions based on intuitive physical reasoning (explained in 
the Introduction) and on the numerical study of (l6| . 

The conclusion is that, at least on a very rough quali- 
tative level, it is sensible to consider R/L^ to be a func- 
tion primarily of q/e and te. Since the q profile tends 
to change more slowly in tokamaks than other equilib- 
rium profiles [58j |. it may be useful to think of a critical 
curve R/ Lx^djE) [111 parametrized by q/e the lat- 
ter quantity containing the essential information about 
the nature of the magnetic cage confining the plasma. 

Collisionality dependence. Even though the de- 
pendence of R/L^ is not as strong as q/e and 7b, 
some discernible inverse correlation between R/L^ and 
v*i might be expected because zonal flows, believed to 
suppress turbulence [60| . should be more strongly 
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FIG. 1. The dependence of R/Lt { (color) on q/e and 7e = 
nrU'^/vthi, showing (a) the individual data points and (b) 
mean values of R/L^ within rectangular bins. 



damped at higher ion collisionality [61 - 64 1 . To isolate 
this dependence, we selected data points for approxi- 
mately fixed te € [0.7,0.8] and q/e € [5,6] (the largest 
number of data points could be found within these nar- 
row ranges and no measurable correlation between R/L^ 
and q/e or te was present). The resulting Fig. ^confirms 
a degree of inverse correlation between R/L^ and v*i. 

Turbulent heat flux. The turbulent ion heat flux 
through a given flux surface is (very approximately!) 
nTiXi/Lxi, where the effective turbulent diffusivity 
c/B)ip/£ y is the (radial) fluctuating 
E x B velocity, r c its correlation time, £ y its poloidal 
correlation scale and ip the fluctuating electrostatic po- 
tential. The latter can be estimated from density fluctu- 
ations using the approximation of Boltzmann electrons: 
e(p/T e w 5n/n (e is the proton charge, n and Sn the 
mean and fluctuating density, respectively). Both the- 
ory of ITG turbulence 12| and the BES measurements 



Xi ~ Su 2 t c , Su 



in MAST [lfj suggest that r c ~ T +i = lyL^ / 'vthipi (the 
drift time; pi is the ion Larmor radius). Collecting all this 
together, we estimate the gyro-Bohm-normalized turbu- 
lent ion heat flux: 



Qj,turb 2i.(R\ fll—\ 



QgB 



Pi 



Ti n 



(2) 



where Q gB = nTiV thi p 2 / R 2 . 

Since ion-scale density fluctuations in MAST can be 
measured directly by the BES system, Qturb can be 
obtained independently of any transport reconstruction 



models such as TRANSP |65j. The method of deter- 
mining Sn/n and £ y using the BES system on MAST 
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FIG. 2. R/L Ti vs. v*i for fixed j E G [0.7, 0.8] and q/e £ [5, 6]. 



(8 radial x 4 vertical channels with spatial resolution of 
2 cm jiH) is explained in detail in [h|. This is done 
from the covariance and correlation functions of the pho- 
ton intensity fluctuations, averaged over the same 5 ms 
intervals for the same 39 discharges as the equilibrium 
quantities studied above, although not in all intervals 
there was good BES data. Restricted to the radial range 
0.6 < r/a < 0.7, the number of available data points for 
Sn/n and l y was 102. 

Inverse correlation between R/L^ and Qturb- It is 
shown in Fig. 02(a) (Qturb vs. q/e and 7 B ; cf. Fig. [Q) 
and Fig. E](b) (Qturb vs. R/L^) that smaller Qturb is ob- 
served where R/L^ is large and vice versa. This is per- 
haps counterintuitive as one might expect that the larger 
R/LTi, the more turbulent the plasma and so the larger 
the turbulent heat flux. This would indeed have been the 
case had R/L^ been externally fixed, as in local flux- 
tube simulations, where Qturb does increase with R/L^ 
121 ljj ]. In contrast, in the real plasma, a certain amount 
of power flows through a flux surface and, if transport is 
stiff, the temperature gradient (along with other equi- 
librium quantities) adjusts to stay close to the critical 
gradient defined by the manifold R/ LT ijC {lE,q/ s). In- 
deed, in plasmas with both power and momentum injec- 
tion, there is a regime with R/L^ close to R/Lt^c where 
the turbulent and the neoclassical (collisional) transport 
are comparable and in which a larger heat flux results 
in lower R/Lt; ■ For a detailed explanation, we refer the 
reader to [59( (see also Fig. 4(b) of 24 1, which should be 
compared to our Fig. |3tb)). In brief, in the neoclassical 
regime, the momentum transport is much less efficient 
than the heat transport, while in the turbulent regime, 
they are comparable (the turbulent Prandtl number is 
order unity HI [H, [37 



66j while the neoclassical one is 
small |28|), so, as a larger heat flux takes the system 
(slightly) farther from the marginal state, this leads to 
much more efficient momentum transport, hence smaller 
velocity shear je , hence a regime with less suppression of 
turbulence and smaller R/Lt^c (see Fig.[3][c), where the 
correspondence between larger Qturb, lower R/L^ and 
lower is shown at approximately fixed q/e G [9, 11]; 
cf. Fig. 3(b) of f59|). Wc stress that all of this happens 
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FIG. 3. (a) Qturb (color) calculated from BES data according 
to Eq. © vs^r/e and 7s (cf. Fig.[Lj); (b) Q tU rb vs. R/L n (d. 
Fig. 4(b) of [2J]); (c) Qturb (color) vs. R/L^ and for a fixed 
range of q/e £ [9, 11], indicated by the dotted horizontal lines 
in (a); open squares are data from Fig. [TJa) for which BES 
measurements were not available (cf. Fig. 3(b) of [591] 'l . 



quite close to marginality and so our experimental obser- 
vation of an inverse correlation between i?/Ly i and Qturb 
provides circumstantial evidence that R/L^ in MAST is 
indeed close to its critical value [dT} . 

Conclusion. Wc have found that the normalized in- 
verse ion-tcmperature-gradient scale length R/L^ has 
its strongest local correlation with q/e and the shear in 
the equilibrium toroidal flow: R/L^ increases with in- 
creasing shear, which is a well known effect, and with 
decreasing q/e, which corresponds to an increasing ratio 
of the perpendicular to the parallel shearing rates. We 
note that a similar dependence of R/L^ (q/e, je) is also 
observed in JET [(58|, suggesting that the inverse corre- 
lation between R/ Lx t and q/e is perhaps ubiquitous, as 
would be the case if R/L^ were gen erally fixed at some 
locally determined critical value [16j. Furthermore, we 
have found an inverse correlation between R/L^ and the 
gyro-Bohm-normalized turbulent heat flux (estimated 
via direct measurements of density fluctuations) and ar- 
gued that this is consistent with R/ always remain- 
ing close to a critical manifold R/L^.c (<?/£, 7sL separat- 
ing the turbulent and non-turbulent regimes [24I [25L (59[ 
(stiff transport). It is thus plausible that we have essen- 
tially produced this critical manifold for the MAST dis- 
charges we investigated. Practically, our results suggest 
that R/L^ can be increased by lowering the q/e, which 
is relatively easier and less expensive than increasing the 
shearing rate in tokamak operations [r39| . 
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